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Eigenvalues of the hermitian Wilson-Dirac operator and chiral properties 
of the domain- wall fermion* 
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Chiral properties of QCD formulated with the domain-wall fermion (DWQCD) are studied using the anomalous 
quark mass m 5q and the spectrum of the 4-dimensional Wilson-Dirac operator. Numerical simulations are made 
with the standard plaquette gauge action and a renormalization-group improved gauge action. Results are reported 
on the density of zero eigenvalue obtained with the accumulation method, and a comparison is made with the 
results for ms q . 



1. Introduction 

Formulation of chiral fermions on the lattice 
has been one of long-standing problems in lattice 
field theories. Several years ago, the domain-wall 
fermion (DWF) formalism W^, which is a Wil- 
son fermion in D + 1 dimensions with Dirichlet 
boundary condition in the extra dimension, has 
been proposed as a new formulation of lattice chi- 
ral fermion. In the limit of large extra dimension 
size, N s — > oo, the spectrum of free domain- wall 
fermion contains massless modes at the edges in 
the extra dimension. 

While the massless modes are shown to be sta- 
ble in perturbation theoryfjUl, their existence 
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may be spoiled non-perturbatively in the pres- 
ence of dynamical gauge fields. We studied this 
issue through an anomalous quark mass m§ q in 
Ref. This quantity measures the magnitude 
of chiral symmetry breaking with the domain- wall 
QCD (DWQCD). 

In this article we make a status report of our 
attempt to understand the results on the N s - 
dependence of m^ q obtained in Ref. || through 
measurements of the eigenvalue distribution of 
the 4-dimensional Wilson-Dirac operator. 

2. Chiral Properties of DWQCD 

We define the anomalous quark mass 

r Ex(^x)P b (o,o)) m 
mB ' = ^Ex(^(*,x)i*(o,o)>- (1) 

This quantity measures the chiral symmetry 
breaking effect in the axial Ward-Takahashi iden- 
tity: 

(V»A°(x)P»(Q)) = 2m f (P a (x)P b (0)) 
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Figure 1. m^ q as a function of N s at a" 1 ~ 1 
GeV and too = 1.7 and at a -1 ~ 2 GeV and 
too = 1.8, where Too is the domain- wall height, for 
the plaquette (f3 = 5.65, 6.0) and RG-improved 
(J3 = 2.2,2.6) gauge action f|. 

+2(J^(x)P b (Q))+i(S a x P b (0)) , (2) 

where A"(x) is the axial- vector current, P a (x) 
is the pseudoscalar density, and J$ q {x) repre- 
sents the explicit breaking of chiral symme- 
try. On a smooth gauge field background, the 
anomalous contribution (j£ q {x)P b (y)) vanishes 
as exp(— cN s ) at large N s ||. Therefore m 5q also 
vanishes exponentially in this case. 

In Ref. H , we carried out quenched simulations 
to study the N s dependence in detail. Two val- 
ues of lattice spacing, a -1 ~ 1 and 2 GeV, are 
explored, using both the plaquette and an RG- 
improved gauge actions. 

Our main results for m§ q are summarized in 
Fig. [Ij The results for the plaquette action arc 
obtained at = 5.65 (a" 1 ~ 1 GeV) and 6.0 
(2 GeV), and those for the RG-improved action 
at (3 = 2.2 (1 GeV) and 2.6 (2 GeV). Solid lines 
are fits to c + ae~^ Ns , and dashed lines to ae~^ Ns 

Important points to note in Fig. [l] are: (i) com- 
paring results for the spatial sizes 12 3 and 16 3 
at a -1 ~ 1 GeV, we find the finite volume ef- 
fects in ui§ q to be small, (ii) tos 9 decreases with 
the lattice spacing, (iii) the magnitude of m$ q 
is smaller for the RG-improved action than for 
the plaquette action, and (iv) from data in the 



range of N s we explore, m^q seems to remain non- 
zero in the limit N s — > oo, in all cases except 
at j3 = 2.6 for the RG-improved action. If con- 
firmed with studies at larger values of N s , the last 
point means that DWQCD realizes chiral sym- 
metry at a -1 ~ 2 GeV only for the case of the 
RG-improved action. 

3. Eigenvalues of the hermitian Wilson- 
Dirac operator and chiral property 

Chiral symmetry of DWQCD can be studied 
also through the transfer matrix in the direction 
of the extra dimension j7|||. When the transfer 
matrix has a unit eigenvalue, chiral symmetry is 
not realized in DWQCD because the left and right 
chiral modes on the two edges in the extra dimen- 
sion couple with each other. 

A unit eigenvalue of the transfer matrix is in 
one-to-one correspondence with a zero eigenvalue 
of the hermitian Wilson-Dirac operator defined 
by 

H w (m ) = 75-CV(-too) , (3) 

which is much easier to calculate. Here, 
Dw{—mo) is the four dimensional Wilson-Dirac 
kernel with a bare mass —mo. Therefore, a failure 
of exponential decay of m§ q would result if Hw 
develops a zero eigenvalue. 

We calculate eigenvalues of by the Lanc- 
zos method using 50-100 configurations at sev- 
eral values of coupling in the range a -1 ~ 1-2 
GeV using both plaquette and RG-improved ac- 
tions. The results from the Lanczos method are 
checked by the Ritz functional method for Hw 
We also study the dependence on the lattice size. 
The maximum lattice at a -1 ~ 1 GeV is 12 4 for 
both actions, while the one at a -1 ~ 2 GeV is 
24 4 for the RG-improved action and 16 3 x 32 for 
the plaquette action. 

3.1. Eigenvalue distributions 

In Figs. H and 3 we plot Monte Carlo time his- 
tories for the six lowest eigenvalues of for the 
plaquette and RG-improved actions. In each fig- 
ure the left panel shows results for a -1 ~ 1 GeV 
and the right panel for a -1 ~ 2 GeV. The lattice 
size at a' 1 ~ 2 GeV is the same as in the previous 
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Figure 2. Monte Carlo time histories for the 
six lowest eigenvalues of H^y obtained with the 
plaquette gauge action. 

work of m 5q shown in Fig. [l[ Open squares plot 
the minimum eigenvalue A, 2 nin and filled diamonds 
are the five higher eigenvalues. 

There is a clear trend that the minimum eigen- 
values become larger for smaller lattice spac- 
ings. Another interesting point is that the RG- 
improved action gives larger values of A^ in than 
the plaquette action, which indicates that the 
RG-improved action has a better chiral behav- 
ior. These trends are parallel to the features we 
noted for m^ q in Sec. |2|. 

3.2. Spectral density 

The spectral density of H\y is defined by 

^) = ^^X>(A'-A), (4) 

where the summation is over the eigenvalues of 
Hw- We are interested in the density of zero 
eigenvalues, p(0), since we expect this quantity 
to be related to the existence of unit eigenvalue 
of the transfer matrix. To calculate this quantity, 
we adopt the accumulation method proposed in 
H), which is based on the relation 

A 2 

A(X) = [ rfA' 2 p(A' 2 ) = -|— £ 1 (5) 

J ° |A'|<A 



Figure 3. The same as Fig. || obtained with the 
RG-improved gauge action. 




dX'p(X') ~ 2 j o(0)A + 0(A 2 ), (6) 



where p{X 2 ) is the spectral density function for 
H%/- We note that, for the small-A expansion of 
j4(A) in (^J), analyticity of p(X) at the origin is 
assumed. 

In Fig. |], we show typical examples of the ac- 
cumulation A(X) from the eigenvalue distribution 
of for the case of the RG-improved action. 
Results for p(0) obtained by a linear fitting fol- 
lowing (^), normalized by the string tension, are 
summarized in Fig. 

Our results for p(0) for the plaquette action are 
consistent with the previous data by Edwards et 
al. (|]. Results for the RG-improved action show 
a similar (3 dependence. A significant difference 
is that the RG-improved action leads to much 
smaller values of p(0) than the plaquette action, 
roughly by an order of magnitude. 

4. Discussions 

We have applied the accumulation method to 
estimate the spectral density at zero eigenvalue 
of the hermitian Wilson-Dirac operator, p(0). We 
found that this method leads to non-zero values 
of p(Q) at a -1 ~ 1-2 GeV for both the plaquette 
and RG-improved actions. 

At a -1 ~ 1 GeV, the non-zero result for p(0) 
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Figure 4. The accumulation A(X) at /? = 2.2 
(a" 1 ~ 1 GeV) and f3 = 2.6 (a" 1 ~ 2 GeV) from 
the RG-improved action. 



is consistent with the finite m^ q in the large N s 
limit observed in Q with both the plaquette and 
RG-improved actions. At a -1 ~ 2 GeV, while a 
consistency also holds with the plaquette action, 
there is an apparent contradiction for the case of 
the RG-improved action since m^ q seems to decay 
exponentially with N s for this case. 

In the accumulation data shown in Fig. |J we 
observe that results are very noisy at f3 = 2.6 
(2 GeV). Since the fit results for p(0) fluctuates 
with volume, it is difficult to determine the size 
dependence. Therefore simulations with larger 
lattices are needed to check if the slope remains 
non-vanishing toward infinite volume. Another 
point to examine is if the analyticity assumption 



Figure 5. p(0)/cr 3 / 2 as a function of (3 for the 
plaquette and RG-improved gauge action. Data 
obtained on the largest lattices are shown. 

for p(A) at the origin is justified if there is a spec- 
tral gap. Further studies are required to clarify 
these points. 
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